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Abstract: 

In this paper, Landau diamagnetism is revisited in the context of nonextensive 
statistical thermodynamics. Exact as well as approximate expressions for partition 
function, magnetization and susceptibility are obtained by using Hilhorst integral 
transformation which provides exact solutions in the framework of Tsallis statistics. The 
results have been reassessed with the standard results to illustrate the effect of 
nonextensivity. 
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1. Introduction 

It has been understood that extensive (additive) Boltzmann-Gibbs (BG) thermal 
statistics fails to study the nonextensive physical systems where long-range interactions 
or long-range microscopic memory are involved or, the system evolves in a 
(multi)fractal space -time. Thus, standard statistical mechanics, consequently, standard 
thermodynamics is not universal and is suitable for extensive systems. A generalized 
formalism is proposed by C. Tsallis which might overcome some of the difficulties that 
nonextensive physical systems exhibit [1,2,3]. 

This generalization relies on a new entropic form for the entropy that is inspired from 
(multi)fractals 

w 

S q =k B ^— (qGR) 
q-l 

where k B is a positive constant and W is the total number of microscopic accessible 
states of the system (for the q<0 case, those probabilities which are not positive should 
be excluded). This expression recovers the well-known standard Shannon entropy in the 
limit q— >1. 

w 

s i =~k B Y,Pi ln Pi 

i=i 

The entropic index q, related to and determined by the microscopic dynamics, 
characterizes the degree of nonextensivity. 

Recently, generalized statistical mechanics has been successfully applied to 
investigate physical systems which exhibits nonextensive features. Amongst them, there 
are stellar polytrops [3], Levy- like anomalous diffusions [4], two dimensional 
turbulence[5], solar neutrino problem [6] velocity distributions of galaxy clusters[7], 
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Cosmic background radiation[8,9] and correlated themes[10], linear response 
theory[ll], thermalization of electron-phonon system[12] and low dimensional 
dissipative systems[13]. 

It should be remarked that what makes S q favourable is that it has, with regard to 
{pi}, definite concavity property for all values of q[14]. 

In order to have a sensible results which are exact or approximate some techniques of 
calculations have been formulated in Tsallis generalized statistical mechanics formalism 
(TT). These are (1-q) expansion [8], factorization method for quantal distribution 
functions [15], perturbative expansion [16], variational methods [16,17], semiclassical 
approximation [18], Feynman's path integrals method in nonextensive physics [19], and 
Green functions [20] . 

In this paper however, we have revisited the Landau diamagnetism where it is 
investigated by incorporating coherent states and using factorization method within 
TT[21]. Here, our motivation is to redrive the physical quantities such as generalized 
partition function, magnetization and susceptibility exactly as well as approximately by 
incorporating Hilhorst Integral transformation in the context of nonextensive 
thermostatistics, moreover, to reassess the obtained results with the standard ones. We 
believe that our approach will help to understand the difference between a 

perfect conductor and a superconductor in terms of entropic index q. 

2. Energy states of an electron in a uniform magnetic field and diamagnetism 

Diamagnetism arises from the quantization of the induced circular current of charged 
particles in an external magnetic field. According to the Lenz rule, the induced dipole 
moments oppose to the direction of the applied field, thus the diamagnetic materials 
weaken it and displacement of the magnetic field out of a diamagnetic material is 
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observed. In case of ideal diamagnetic material such as superconductors the interior of 
the materials is free of the applied field, this effect is known as Meissner-Ochsenfeld 
effect more commonly referred to as Meissner effect (Berlin, 1933). 

Here, the magnetic behaviour of the electrons has been considered. Permanent 
magnetic dipole moments of electrons and nuclei are ignored and only the influence of 
the external magnetic fields on the electrons is taken into account. 

The Hamiltonian of a nonrelativistic free electron in an external magnetic field is 



H = 



l 

2m 



P+-A 



(1) 



where e is positive. The Schrodinger equation 

H x ¥ = eyf 



(2) 



is invariant under the Coulomb gauge transformation where A • A = . We consider a 
uniform external magnetic field H along the z-axis, and choose the vector potential 



A = -L 
2m 



f 1 1 ^ 
— H v -H^O 

2 y ' 2 x 



(3) 



— > — > 

so that H = H k , thus the Hamiltonian of an electron is written down as 



H = 



2m 



1 V / 



-mOty 



+ 



P., + — mCOx 

y 2 



+p: 



(4) 



eH 

where (0 = is the electron cyclotron frequency m is the mass of electron, and c is 

mc 



the velocity of light. Introducing the operators 

f in ^ 



jz + =P ±iP ± 



{x±iy) 



(5) 



where £ = 



f h ^ 



y m(0 j 



is the classical radius of ground state Landau orbit and noting that 



they obey the commutation relation 

It can be shown that the Hamiltonian becomes: 



(6) 



H = 



> 2 f 



2m 



2m 



+ -Ti(D . 
2 



(7) 



Eq.(7) can be separated into parallel and transverse components: Thus, the solution of 
the Schrodinger Equation of the Hamiltonian which is given by Eq.(2) consists of the 
motion along the z direction of a free particle whose Hamiltonian H z and the motion of 
harmonic oscillator whose Hamiltonian H t and move on the xy plane 



H = H Z +H t 



(8) 



where 



H = 



2m 



(9) 



The energy of free particle is known. On the other hand, choosing the convenient 
basis of set, for the harmonic oscillator the energy eigenvalues are obtained: 



eH O 
n + - 
2 



mc 



n=0,l,2, 



(10) 



The energy levels of Eq. (10) are degenerate. The number of states g which belong to 
a discrete n are calculated [26]. The levels of the field-free case which superpose have 

an energy e = (p x 2 + Py ) which lies between the energy interval 2ji B Hn and 
2fi B H{n + \). Thus 



g = \\dP x dP y dxdy 
h J 



(11) 



2n B Hn<£<2n B H(n + l) 
with 

lm Imc 
The integral over x and y just yields the base area of the container, the momentum 
integrals can be solved by substituting plane polar coordinates: 

h i h 



or 

l-8ijL (13) 
he 



where P„ = jAm fJi^Hn 



Degeneracy is independent of n. It vanishes while H->0 , since we have assumed the 
limiting case of a continuos spectrum in Eq. (12). In our case we take a cylindrical 
sample of length L z , radius R, therefore base area is S=tiR 2 and it is oriented along the 
magnetic field H. 

We consider a system, which is in a thermodynamical equilibrium state, and 
fluctuate in terms of energy with its environment, the generalized partition function for 
the canonical ensemble within the TT formalism is given by [1,2]. 



i 



Z =Tr [l-(l-<z)j3H]w (14) 



where (5 = . The generalized free energy F of the sample is of which the electron 
k B T 

density n ,is obtained from the generalized partition function Z , 



and generalized magnetization M q is 



dF 



dH 



M q =~^± (16) 



and generalized susceptibility is 

*< = -is- (17) 



3. Hilhorst Integral Transformation and Generalized Canonical Partition 
Function For Landau Diamagnetism 

Hilhorst integral transformation [22] and extension for q<l bridges the generalized 
statistical mechanics to its respective standard quantities thus generalized 
thermodynamics could easily be established [23]. Therefore the generalized statistical 
mechanics can be regarded as a Hilhorst integral transformation of respective standart 
quantities. From the definition of Gamma function [24] 

CO 

T(a)= ^dte-'t"' 1 (18) 

o 

substituting 

t = v[l-{l-q)(3H] and a = — intoEq.(18) 

l-q 

one obtains: 

1 1 cc _j 1 

[l-(l-q)(3H]vr q =— -$dvv q ~ cxp[-v{q-lW] (19) 

— 1 o 



On the other hand, the Trace of Eq.(19) provides the generalized Canonical-partition 
function Z ? for Landau Diamagnetism, which can be evaluated by choosing a 

convenient set of basis states: 



1 r 1 



(20) 



where 



(21) 



Substituting the components of the Hamiltonian of Eq.(9) into Eq.(14), the partition 
function Z v (v(q -l)jS) is separated in two parts 

Z 1 (v(q -l)fi)= Z u/ {v{q - l)p )Z l± {v{q - l)fi ) (22) 
where Z y/ is the partition function, which is parallel to the cylinder, 



Z y/ (v(q-l)(3)=L z 



m 



27Ch 2 (3v(q-l) 



(23) 



where L z is the height of the cylinder. On the other hand the transverse part is 

&xp(-v(q-l)x) 



Z l± (v(<l-l)P)=g 



l-exp(-u(g-l)jc) 



(24) 



eH 



where g is the number of the states which is given by Eq.(13) and x = f5h . Using 

mc 

Eq.(13) for the geometry of the cylinder the number of the states is 



g = 2KR' 



eH 

he 



By substituting of Eq.(23) and Eq.(24) into Eq.(22) one obtains 
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' mCO N 


' m ^ 


2ith j 


K 27th 2 , 



fe-i)^r 



— 1 



where 



1 -Ihi 



J [l-exp(-u(q-l)x)] 



(25) 



(26) 



One can write Eq.[26] as 



/„=E fW ,_1 2 ^ exp(-u[(l -(1 -?)*(»+ X)J 



(27) 



In order to terminate the integral transformation, one may, once more, recourse to the 
change of integral variable x> in Eq.(27) 



u =v 



l-(l-q)x 
then Eq.(27) reads, 



n + - 
v 2 



(28) 



/„ =r 



' l P 



q-l 2 



E 

)n=0 



l-(l-q)x 



n + - 

v 2 y 



i i 

— +— 

l-q 2 



(29) 



Since — — > one can draw the conclusion l<q<3 which is in agreement with the 

q-l 2 

distribution function of paramagnetism where a fractal approach is followed [25]. 
Finally, Eq.(25) becomes 



f m(0 Y 



27th 



m 



V~"" 1 A^'" 4 J 



2lth' 



n=0L 



l-(l-q)x 



n + - 



v 2 yj 



i i 

— +— 

l-q 2 



(30) 



where r{a,b)) = ^^e b ^ a 

r{a) 



1 h 1 

a = ,b = — 

q-l 2 
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When a— ►co i.e. q — >1 it is easy to show that Z q (P) turns to the standard partition 
function, i.e. 



f m(0 ^ 



m 



(31) 



since 



limr(a,^)^l [27]. 



In the case of high temperature i.e x = (5ha>«l generalized partition function 
simplifies and Eq.(30) becomes 



2 ^U^ 2 J l-[l-{l- q ) x ]-^\ 



q>l 



(32) 



since l-(l-q)x 



( n 

n + — 

v 2 ; 



exp 



(l-<7> 



f n 



n + - 



including zero-point energy. 

Now, Let us obtain the generalized partition function for the case q<l.The contour 
integral, which is suitable for q<l case is [ 27 ] 



Y{a)—Ut{-t)- a e- t =1 
2k j 



(33a) 



where a>0 



Substitution of 



t = v\l-(l-q)(3H]; and a = l + — *— (a>0 since q<l) (33b) 

l-q 
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into Eq.(33a) leads to: 



[i-(i-^K]w=-^r 

171 



1-9 



f t/v(-u)" 1_ X-? exp(-u)exp(u(l-^) J SK). (34) 



In order to obtain the generalized partition function from Eq.(34) a convenient set of 
basis states are chosen, then the trace of the corresponding matrix is taken. Thus, the 
partition function of the nonextensive system is obtained for q<l: 



z,G8)=r 



2-q_ |j_ 



§dv(-v) \- q 1 exp(-u)z 1 (-jSu(l-?)). 



(35) 



Similar to the arguments which is given above for the q>l case ,the components of 
Zj (- (3l)(l-q)) consists of two parts one is parallel other is transverse to the axis of the 



sample. Thus, 



(36) 



where 



Z v/ {-v{l- q )/3)=L z 



m 



27Zh 2 {-v{\-q)P) 



m 



v (37) 



and 



exp(-v(l-^)^) 
1 - exp(- v(l - q)x) 



(38) 



The number of the degenerate states g in Eq.(38) is given by Eq.(13) . 

By substituting Eq.(37)and Eq.(38) into Eq.(36) the partition function which given 
by Eq.(36) is found. When Z { {-^v{i-q)) given by Eq.(36) is replaced in Eq.(35) 
Z (j8) is obtaied as : 
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zM=v 



f mCO V 



m 



A/2 



27thp 



l-q 



+ 1 



(i-^r 



f 1 3^ 

- + - 



n=0 



l-(l-q)x 



n + — 

v 2 yj 



i i 

— +— 

l-q 2 



^1-9 2 

Finally,in view of Eq.(30) and Eq.(39) Z (j8) is written for different q's: 



m 



J 



1 1 

~q-\ 2 



(g-i)Xr 



,9-1 
+1 



(1-?)/ 



<?J/2r 



1 3 
+ 



l-q 2 



q>\ 



q<l 



(39) 



(40) 



where 



n=0 



!-(!-«?) 



n + - 



i l 

— +— 

1-0 2 



(41) 



It is wise to take the logarithm of Z q (P) and its derivative with respect to x since one 
needs them for advancing further: 

lnZ ? =lnC + ln;t + lnX 9 (42) 



where C stands for all constants and 



— In Z =- + — lnX 

dx x dx 



(43) 



4. The generalized magnetization and susceptibility of Landau Diamagnetism 

Materials are generally classified by the sign and the magnitude of their 
susceptibilities. The majority of the materials (solids) are diamagnetic for which 
susceptibility % is negative. The generalized magnetization M q in terms of the free 



n 
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energy is given by Eq.(16). On the other hand, when the relation between F q and Z q is 
taken into account the generalized magnetization follows: 



M„ =n — Z l ; q —\nZ 
mc dx 



(44) 



since 



dx _ p eh 
dH mc 



(45) 



and 



dx /3 dx 



m Z 1 -^i nZ . 



(46) 



Substituting the derivative of Z q with respect to x which is given by Eq.(43) 
into Eq.(46) one obtains for the generalized magnetization 



mc 



l l 



x 



V ) 



oo 

I 


( n 


l- 


(i- 


q)x 


( n 

v 2 1 


2 


n=0 


v 2 ; 













(47) 



which is is the generalized exact magnetization expression for the diamagnet, since 



3.1 ' 2 ^ 



n=0 



n + 

V 2 A 



l-(l-q)x 



n + — 

v 2 y 



i _i 

l-q 2 



X. 



(48) 



In the high temperature limit, it can be shown that Eq.(41) and Eq.(42) could may 
take approximated values, i.e. 



' ~l-Y 2 

1 



(49) 



and 
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— lnX„ = 

dx q 



l-y„ 2 



(50) 



where 



1 i 

— +— 

l-q 2 



q 2 ; Y q 2 =[l-(l-q)x]^ q + 2 , Q q =^-\nY q . (51) 



Thus, the approximated value of magnetization which given by Eq.(47) is obtained in 
terms of Y q 

f 1 , v2\ 



M„ = 



nfie 



mc 



1 



- + 



1+5 

l-ri 



z i- 9 



(52) 



nfie 



mc 



x 



1 1 



Yi 



A 



Z 1 " 



i-(i-<zV 



(53) 



since 



Qa=~ 



1 1 Yl-fl 
-+ 1 



[l-q 2 



1 



(54) 



The magnetization is found to be a standard one when q— ►! i.e 



M 1 = 



nfie 



mc 



11 , x 

coth — 

x 2 2 



(55) 



since 



1 

2 



yj 9=1 = exp(-^); Z^U^l 



(56) 



, x e 1 +e 2 
coth — = 

2 * 

e 2 — e 2 



(57) 
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We keep only the lowest-order contribution of x for Mi in Eq.(55) then the usual 
expression is obtained 



nhe 
Ylmc 



-x , 



(58) 



Now let us consider the diamagnetic susceptibility: 

In accordance with Eq.(17) the magnetic susceptibility per electron is given by 



1 _ he dM 

X q = — P 



n mc 



dx 



(59) 



dx eh 

since = /3 — in accordance with Eq. (45). Substitution Eq.(44) for M q into Eq.(59) 

dH mc 



'he ' 












_x + dx X «_ 


^ mc j 


3bc| 





leads to 

x q =-P 

where Eq.(43) is taken into account. Rearranging, then it becomes 

Xa=~P 



(60) 



'he} 2 
y mc j 



dx z 



dx 



dx 



(61) 

Under the approximation, which is made for diamagnetization, and along the line 
that M q is derived, susceptibility is obtained in terms of Y q 



he 



y mC j 



(q-2) — + 



l 4Y: 



Q + 



x 2 ' 1-F 2 q 



l + Y. 



2 V 



l-Y: 



2(«-l) 



X 



i-y: 



2 \ 



2 A 



i-y: 



Q 2 

&q 



A-q 



(62) 



It is easy to verify that in the q— »1 limit, Eq.(62) leads to the standard diamagnetic 
susceptibility. 
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Xi=-P 



he 



y niC J 



- + 



(63) 



or 



Xi=~P 



f he A 



y MC J 



1 1 ]2 x 

— + — esc h — 
x 2 4 2 



(64) 



since 



Y{=e- X \ 



csch — = 
2 



(65) 



As in the magnetisation one may keep the lowest order contribution in x for 
susceptibility then Eq.(58 ) becomes the well-known standard expression [28]: 

s2 



P 



f eh ^ 



X^-^r — ■ (66) 
3 ^ 2mc 

Conclusions 

As is stressed in the introduction BG statistical mechanics is based on the extensive 
(additive) entropy that optimisation of entropy under the constraints yields an 
exponential dependence on energy. We study Landau Diamagnetism within the context 
of generalised TT in which long-range interactions, long range memory multifractal 
space-time and so forth, i.e. nonextensivity is taken into considerations. In TT 
formalism entropic index q play a major role and it is an indicator of nonextensivity. 

In an earlier paper [21] Landau diamagnetism was considered in the limit where 
factorisation approximation [15] is valid. Here, we employed the Hillhorst integral 
transformation in order to obtain the generalised partition functions Z q exactly,then 
magnetisation free energy F q , generalised magnetization, M q and finally generalized 
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susceptibility J£ are calculated, high temperature limits of those quantities are also 
discussed. 

In the q— >4 limit, standard usual expressions for Zi,Fi,Mi,%j are also indicated. 

According to our justification this paper is on the line to obtain the exact solutions 
for the nonextensive thermostatistical systems, one of them is Landau diamagnetism. 
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